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A new LES model derived from 
generalized Navier-Stokes equations 
with nonlinear viscosity 
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Abstract 

Large Eddy Simulation (LES) is a very useful tool when simulating turbulent 
flows if we are only interested in its “larger” scales (see [II], (9l or O, for ex¬ 
ample). One of the possible ways to derive the LES equations is to apply a filter 
operator to the Navier-Stokes equations, obtaining a new equation governing the 
behavior of the filtered velocity. This approach introduces in the equations the 
so called subgrid-scale tensor, that must be expressed in terms of the filtered ve¬ 
locity to close the problem. One of the most popular models is that proposed by 
Smagorinsky, where the subgrid-scale tensor is modeled by introducing an eddy 
viscosity. 

In this work, we shall propose a new approximation to this problem by applying 
the filter, not to the Navier-Stokes equations, but to a generalized version of them 
with nonlinear viscosity. That is, we shall introduce a nonlinear viscosity, not as 
a procedure to close the subgrid-scale tensor, but as part of the model itself (see 
below). Consequently, we shall need a different method to close the subgrid-scale 
tensor, and we shall use the Clark approximation, where the Taylor expansion of 
the subgrid-scale tensor is computed (see O and QSl). 

Keywords: Large Eddy Simulation, nonlinear viscosity, Clark approximation. 


1 Introduction 

It is generally accepted that Navier-Stokes equations accurately model the behavior 
of incompressible viscous fluids. Simulating numerically these equations is possible 
today thanks to high computing power available. However, Direct Numerical Sim¬ 
ulation (DNS) is limited to relatively low Reynolds numbers, because to simulate 
Navier-Stokes equations for a given value of Reynolds number (Re) are needed, at 
least, 0{Rt^^) degrees of freedom (see ||9ll). 
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There are several methods to overcome this problem. For example, Reynolds Av¬ 
eraged Navier-Stokes models, k — e models or Large Eddy Simulation (LES). In this 
paper we shall talk about LES. Usually, LES is presented as an “averaged” or “filtered” 
version of Navier-Stokes equations. 

Let us consider a “filter” operator / —>^ / (a spatial filter, a time filter or both), and 
let us assume that it is linear and commutes with spatial and time derivatives (see HI 
for different examples of this kind of filters). If we apply this filter to Navier-Stokes 
equations 

1 

^-l-(Vu)u-|-Vp-:/Au = f, (1) 

at po 

V • u = 0, (2) 


(where u is the velocity field, p is the pressure, and f is the density of external forces), 
we obtain 


du , _ 1 

— —I- (Vu) u H-Vp — i/Au = f — V • T, 

at Po 

V-u = 0, 

Tj^j j Uj Uj^ llj , 


(3) 

(4) 

(5) 


where r is the so called subgrid-scale tensor. To close the LES model is necessary 
to express the subgrid-scale tensor only in terms of u. One of the most popular LES 
models is that proposed by Smagorinsky, where the subgrid-scale tensor is modeled by 
introducing an eddy viscosity Vf. such that 


T = -2l/eD, Ve = 

iCsAxfV2\t)l 

(6) 

D = 1 (Vu -f Vu^) , 

/ 3 _ _ 


|D| = 1 1 , 

(7) 


Ax is the subgrid-scale characteristic length and Cs is a constant chosen to allow the 
model to emulate the kinetic-energy dissipation predicted by Kolmogorov (see El, El 
orQOl). 

The Smagorinsky model is sometimes too dissipative, and it has other theoretical 
problems, as the fact that subgrid-scale tensor in is odd in u, when by definition 
(see 0), it must be even in u. This happens because 0 does not try to approximate 
the subgrid-scale tensor as defined in Q, but (as we have said before) to emulate the 
kinetic-energy dissipation predicted by Kolmogorov. 

Interested reader can see other LES models in El, 03, HI, na or El- 


2 Deriving the new LES model 

In the previous section, we have pointed out that the Smagorinsky model introduces 
dissipation through (j^, in a way that these equations are more a model than an approx¬ 
imation of the subgrid-scale tensor. 


2 



In this section, we shall propose to model this dissipation by introducing a nonlinear 
viscosity in the Navier-Stokes equations and then we shall apply the filter. Thus acting 
we introduce the dissipation predicted by Kolmogorov directly in our model and, when 
applying the filter, we shall close the LES model by approximating the subgrid-scale 
tensor using Clark approximation (see subsection |2.2| i. 

2.1 Introducing the nonlinear viscosity 

Let us consider the following generalization of Navier-Stokes equations 

Po + (Vu) = pof -p V • T, 

V • u = 0, 

where the stress tensor T is given by 

T = -pI-p2^eD, D = i (Vu-P Vu^), (10) 

and where the dynamic viscosity po has been substituted by the effective viscosity ptg 
depending on the norm of the strain rate tensor D. To fix ideas, let us choose the 
following effective viscosity 

/re = Me(|D|)=/xo(l + A"|D|2)« (11) 

where g > 0 and A > 0. 

Remark 1 If we take g = 0 or A = 0, equations ([8|)-([TT]) are the Navier-Stokes 
equations. Otherwise, they are known as Ladyzhenskaya model, and a unique global 
weak solution is guaranteed if g > 1/4 (see j^)- We obtain the Smagorinsky model if 
we take q = 1/2, and properly choosing A, we can recover the dissipation predicted by 
Kolmogorov. 

Lrom the above remark, let us consider the effective viscosity given by GD with 
g = 1/2 in what follows. Thus, we introduce the dissipation predicted by Kolmogorov 
as part of the model, and we shall approach the subgrid-scale tensor in another way 
(see @-@). 


( 8 ) 

(9) 


2.2 Filtering the equations 

In this subsection, we apply a filter to the above equations (([8|)-([TT])). We are interested 
in linear filters that commute with spatial and time derivatives. Although there are 
several possible choices (see M), we shall consider here a Gaussian filter; 





t,y- x)/(s,y) dy ds 


with 


G(s,y) 


1/2 3/2 

iL 


/ 7TS^+7L|yP \ 

V w ) 


( 12 ) 


(13) 
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where 77 > 0 is a small parameter related with the size of the filter, and 7 ^ > 0 , 7 l > 0 
are parameters related with the shape of the filter. 

Following the ideas of Clark (see 13 and 03), we approximate / in ( [T^ by its 
Taylor series, and then we can easily prove that this filter verifies for small values of p: 


f = f + V^ —S + —A/ 

_7t ot‘^ 7l 


f = f-v^ 




_7t dt^ 7i 

fg = f9 + 2 ??^ 




1 dfdg 1 ^ J 

— - V/-Vp 

7 t ut at 


07 


f9 = f9 + 


+ f 


2 dfdg 1 ^ 

— - 

7t dt dt 7l 
1 d'^g 1 


(14) 

(15) 

(16) 
(17) 


^7t dt^ 7l 

and for any smooth enough real function /?, 


Ag 




l3{f)9 = P{f)9 + vHW{f) 


1 dfdg 1 - 

— - 

7 t or or 71 , 


(18) 


+ /?"(/) 


1 


7t V 


1 


+ —V/-V/ 

IL 


0 ( 7 " 


If we apply filter (12 1 -( 13 1 to equations ([^-( 11 1 (with q = 1/2), and we take into 
account ([T4|)-([r8|), we obtain 


O — -I 

+ (Vu) u H-Vp = f + V • (S - r). 

dt pq 

V-u = 0, 


where t is given by 


— Uj/LLg UfU 


= 2 


2 Ujj Uji LLj 

1 Oui duj 1 „ „ 

-+ — Vu* • Vu^- 

_7t at 7 l 


rf + OCt?"*), 


and S is given by 


= - 


2 

Po 
= 2v 


{pe{mT> 

(1 + + A 71 + 

D,i 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 


K 2 Y 


- 1/2 r> „2 


--(l + A^Ki 


+ 2A7l + A"iTi) ^ + 
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where we have introduced the notation 


V = Ki = |D 

m,n—l 


1 

7t 


dDr, 


dt 


H- 

IL 


VDr, 


K. = — 


1 

IT V 
1 dKidDij 


7t dt dt 


—VitTi • VKi, 
IL 

i-V-ftTi ■ V-Di, 

IL 


(24) 

(25) 

(26) 
(27) 


Remark 2 Formulas ([23|-p7]l remember us those obtained in the dynamic procedure 
proposed by Germano (see Q and mi), so this work could be seen as a new way to 
deduce a similar kind of models. 


3 An example: application to Burgers equation 


In order to evaluate this model in the simplest possible case, we are going to apply, in 
this section, the same method to Burgers equation. The Burgers equation 


du du d'^u . 


(28) 


is often considered as an one-dimensional version of the Navier-Stokes equations. 
Next, we shall introduce a generalization of Burgers equation 


Po 


du du 


Pof + Po 


A 

dx 



(29) 


that can also be considered as an one-dimensional version of 


3.1 Filtering Burgers equation 

Let us now define the equivalent filter to ([T^-([T3]l, but one dimensional in space, so 
we can apply it to the generalization of Burgers equation (291, 


fit,x) = 


n 

—oo j —oo 


G(s -t,y- x)f{s, y) dy ds 


(30) 


with 


1/2 1/2 


7ts + ILV^ 

47^2 


(31) 


If we now apply filter ((30ll-((3T]i to equation ( |29l l, and we take into account the 
formulas corresponding to ([g)-0 in the one-dimensional spatial case, we obtain 
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( 32 ) 


du _du - d 


where 


1 


^ 2 


r_ n 


’ 1 

fdu^ 

1 

^ du\^ 

u^ — (w)^ 


7t 

My 

^Tl 

\dx) 


V'^+Oir,^), (33) 


^ 9\ 1/2 

duV\ du 


® = + I i=-'[^* + A=AV]|+0(,‘), (34) 


and 7^4, are defined by 

7^4 = ( 1 + 


^5 = 1 1 + AM — 
ox 


1 ^ 


V 2\ 1/2 

dx ^ 

V 9\ -3/2 
du^ 


j_ /^_^V 

7t V dtdx ) 


(35) 

(36) 


7l 


3 + 2AM — 
dx 


du'' 


3.2 Numerical comparison with some analytical solutions. 

Let us consider the following family of analytical solutions of the Burgers equation 

—2v d(j) 


u = 


(j) dx 


(37) 


where 


(j) = Aq + Bqx + [Ai sin((^ix) + Bi cos(iu;ia;)] exp(—z/wjf) ( 38 ) 
+ [A2 sin(a;2a;) + B2 cos(a;2a:)] exp(— 

Next, we shall use MacCormack scheme (see ifT^ l to solve numerically problems 
and ([3^-(|36|l and to compare the obtained results with the analytical solution (j37]i. 
We shall compare with three different analytical solutions, obtained choosing the 


following three different sets of values for the constants of (37 1 -( 38 1 , p and A, 


wi = 27r,W2 = 507r,^o = 1000, Bq = -10, = 7,^2 = tA-, 

TT IOOtt 

Bi= B2=0 ,v = ^ = 0.1, A = 1.8e6, 

3 3 

U)i = 2 tt,uj2 = 507 r,>lo = 10, Sq = -7,44i = -, ^2 = 77—, 

TT IOOtt 


(39) 


( 40 ) 


Bi = B 2 = 0,1^ = 


1 


50000 


ri = 0.1, A = 5000, 
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( 41 ) 


— 15 1 

OJi = 27r,U!2 = IOOtT, ^0 = 1; = “10^ ; ^2 = : 

TT IUtT 

= —, 132 — -, — -■ 

TT ’ IOOtt’ 50000 


jy = 0.01, A = 500, 


where i/ can be understood here as the viscosity, or as the inverse of Reynolds number 
in the non-dimensional version of ( |28] l. 

We can see in figures [T][^the good behavior of the numerical approximation of (32 1 - 
( [3^ if compared with the numerical approximation of ( [28] l. For these simulations, we 
have chosen in (j3T| Jt = Jl = 6. We have used a spatial step of Ax = 0.1 and a time 
step of At = 0.01. We have plotted in blue the exact solution, in green the numerical 
approximation of the Burgers equation ( [28) l, and in red the numerical approximation of 
(32 1-([3^. We have represented the solutions at f = 1. 
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X 10'® 



Figure 2; Constants (401, Ax = 0.1, At = 0.01, t = 1. 


X 10'“ 



Figure 3; Constants (41 1 , Ax = 0.1, At = 0.01, t = 1. 


Despite we have taken a relatively large discretization step, model (32 1 -( 36 1 gives 
figures [T]|^ a good approximation of the analytical solution of (|28|l, that is, we are 
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able to approximate the “large scale” behavior of the solution without needing a small 
discretization step. 

Next, let us show what happens when we decrease the discretization step. In figures 
and we plot the analytical solution and their approximations for the same set of 
values of the constants as in hgure|^ but with smaller discretization steps (Ax = 0.05, 
At = 0.005 in hgurej^and Aa; = 0.01, At = 0.001 in hgure|^. 



Figure 4: Constants (41 1 , Aa; = 0.05, At = 0.005,1 = 1. 


We can see that the MacCormack scheme has a bad numerical behavior when ap¬ 
proximating Burgers equation (28 1 (oscillations increase when decreasing discretiza¬ 
tion step), while its behavior when approximating equations (|^-(|^ is very good. 

Finally, we plot in hgurej^the same analytical solution (set of values of the con¬ 
stants ( pT] )) and several approximations of equations ([3^-((36ll with different values of 
A and discretization steps. This figure shows how the accuracy of the approximations 
depends not only on the discretization step, but also on the election of A. 
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X 10’’ 



Figure 5: Constants (41 1 , Aa; = 0.01, At = 0.001, t = 1. 


xIO"* 



Figure 6: Constants ([4l]i, different values of A. 
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3.3 Another numerical example 


Now, we shall use again the MacCormack scheme to solve numerically equation (28 1 
with periodic boundary conditions at a; = 0 and a; = 1, a constant initial velocity and 
a periodic external force. We show in figures [7]|^the results obtained taking u(0, x) = 
2.3, v = 1/5000 and the external force / given by 


/ = 2.3 + sin(47ra;) + 2.9 sin(997ra;) + sin(47rf) + 2.9sin(997rf) (42) 

In figure [7] we have plotted different approximations of the solution at point x = 
0.5. We have used four different different discretization steps (Aa; = 0.1, Aa: = 0.05, 
Ax = 0.025, Ax = 10“"^). In order to compare the accuracy of the filtered model, 
let us consider that the “exact” solution is the one given by the smallest discretization 
step (we have plotted this solution in black in figures |7][^. We can see then that the 
filtered solution (computed with the larger step size. Ax = 0.1, and plotted with a dis¬ 
continuous line) is more accurate than the approximations obtained by solving Burgers 
equation (|28|) with smaller step sizes. 



Figure 7: Different approximations at x = 0.5. 
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We have plotted in figure some approximations of the solution, at the instant 
t = 0.7, for the discretization steps mentioned above (Aa; = 0.1, Ax = 0.05, Ax = 
0.025). We can observe that filtered approximations (discontinuous lines) are closer 
to the “exact” solution (Ax = 10“^) than the non filtered approximations (continuous 
lines) with the same, and even smaller, discretization step. This is specially clear for 
the coarser discretization step (Ax = 0.1), where the filtered approximation essentially 
captures the exact solution, while the unfiltered approximation remains far from it. 



Figure 8: Different approximations at f = 0.7. 


4 Conclusions 

A new LES model has been deduced. We have applied a filter to the generalized Navier- 
Stokes equations with a nonlinear effective viscosity and then we have approximated 
the subgrid-scale tensor using the Clark approximation. The new model thus obtained 
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reminds us the dynamic procedure of Germano. 

We have performed some numerical simulations to test this new model. We have 
solved Burgers equation and the model ([32|)-([3^, obtained with the same technique 
used previously for Navier-Stokes equations. We have compared the results with an¬ 
alytical and numerical solutions of the Burgers equations, and we have seen the good 
numerical behavior of the model that we have deduced. Even for large discretization 
steps the hltered model provides quite good approximations of the exact solutions. 
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